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A WIND-TUNNEL EVALUATION OF ANALYTICAL TECHNIQUES 
FOR PREDICTING STATIC STABILITY AND CONTROL CHARACTERISTICS 

OF FLEXIBLE AIRCRAFT 

By Irving Abel 
Langley Research Center 

SUMMARY 

An experimental evaluation of analytical techniques for predicting certain stability 
and control characteristics of a large flexible aircraft is presented. Analytical methods 
based on both the modal approach and flexibility influence coefficients are developed to 
predict the aerodynamic characteristics of a flexible airplane. These methods are then 
applied to a flexibly scaled model of a supersonic transport configuration. Comparisons 
of wind-tunnel data, calculations based on the modal approach, and flexibility influence 
coefficients are presented over the Mach number range from 0.6 to 2.7. An examination 
of the results obtained from this study indicates that both analytical techniques predict 
reasonably well the effect of flexibility on the basic longitudinal characteristics and that 
both techniques give generally comparable results. 

INTRODUCTION 

Effects of aeroelasticity on aircraft design have been considered for many years. 
The aeroelastician has been primarily concerned with such problems as flutter, diver- 
gence, response to turbulence, and control effectiveness. The effect of flexibility on air- 
plane stability has usually been relatively small and could be taken care of by empirical 
corrections to measured or calculated stability derivatives. The effect of aeroelasticity 
on the static aerodynamic characteristics, however, has become a major concern for 
large flexible aircraft that operate in both the transonic and supersonic flight regimes. 

No thorough experimental evaluation of analytical approaches for predicting aero- 
elastic effects or stability characteristics is available. Recent papers by Roskam, 
Holgate, and Shimizu (ref. 1) and by Chevalier, Dornfeld, and Schwanz (ref. 2) present 
some comparisons between theory and experiment, but these are limited by the flexibility 
of the model investigated. Langley Research Center, therefore, in cooperation with The 
Boeing Company, undertook a program to provide a comparison of analytically deter- 
mined and wind -tunnel -measured- rigid and flexible aerodynamic characteristics of a 



proposed supersonic transport configuration. Two independent analytical approaches 
were used: a modal technique by NASA and a direct influence -coefficient technique by 
The Boeing Company. (The analytical approaches are described in appendixes A and B.) 
The experimental portion of the investigation was conducted over a Mach number range 
from 0.6 to 2.7 and utilized separate rigid and flexibly scaled models. 

This paper provides a description of the wind-tunnel models (including the struc- 
tural representations used for the flexible model), a description of the analytical tech- 
niques being evaluated, and a comparison of calculated and measured results. A brief 
overall view of this program is given in reference 3. 

The experimental results presented are not directly applicable to a full-size air- 
plane since they do not include inertia effects associated with deformation of the struc- 
ture due to load factor. The results are, in effect, "massless airplane characteristics" 
since the model experiences a constant (n = 1) load factor and the ratio of gravitational 
forces to aerodynamic forces (Froude number) was not scaled. Both analytical 
approaches, however, can include these inertia loadings for the free-flight case. 

SYMBOLS 

The moment reference center is located at 45 percent of c. 


a speed of sound 






b 


aerodynamic influence -coefficient matrix 

matrix defining deflection on right wing due to a unit load on right wing 

element of aerodynamic matrix [a] which defines pressure coefficient 
at panel i due to a unit angle of attack at panel j 

wing span 



matrix defining deflection on right wing due to a unit load on left wing 


c local chord length measured streamwise 


c reference chord 


c L 


lift coefficient, 


Lift 

qS 


2 



lift coefficient at a = 0° 


lift coefficient due to pitch rate, - ■ , 

d(BL 

\2V, 


lift -curve slope, 


3a lor=0° 


elevator effectiveness in lift, 


rolling -moment coefficient, 


Rolling moment 


effective dihedral parameter, 


aileron effectiveness derivative, 


pitching -moment coefficient, Pitching moment 

qSc 

pitching -moment coefficient at a = 0° 


pitching moment due to pitch rate, 


8&L 

\2V/ 


static stability derivative, 


elevator effectiveness in pitch, 


normal -force coefficient, 


Normal force 


normal -force curve slope, 


9 « Uo° 


yawing -moment coefficient. 


Yawing moment 


directional stability parameter, — — 


side -force coefficient, Side force 

qS 
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C Y/3 

F 

{f} 

g 

hi(x,y) 

H(x,y) 

M 

I 

Ki 

( L ) 

L 

Li 

m 

m (x,y) 

Am 

M 

Mi 

{m} 

Mi 
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side -force parameter, ^ 
air load (see table II) 

matrix defining total resultant force on each wing panel 

acceleration due to gravity 

shape of ith vibration mode of model structure 

vertical displacement of point (x,y) 

indices 

mass moment of inertia in pitch 
generalized stiffness of ith vibration mode 
matrix defining lift on each wing panel 
sum of 

element of matrix {l} defining lift on ith panel 

total number of modes including rigid-body pitch and plunge 

mass distribution of flexible model 

known incremental mass 

airplane total mass; Mach number 

generalized mass of ith vibration mode 

matrix defining pitching moment about reference axis for each wing panel 
pitching moment per unit generalized coordinate 



load factor, ^ lUX 
’ Weight 

normal force per unit generalized coordinate 
total lifting pressure over wing 

pressure distribution due to downwash of ith mode per unit generalized 
coordinate 

1 2 

free -stream dynamic pressure, -pV 
generalized coordinate of ith vibration mode 
generalized coordinate per unit angle of attack, 
ith mode generalized force 


reference wing area 


matrix defining area of each wing panel 


free-stream velocity 


w(x,y) 


downwash on wing 


aircraft total weight 


matrix defining weight of each wing panel 


x,y,z 


distances along coordinate axes 


X,Y,Z orthogonal coordinate system, body axes 


matrix defining distance from each panel to reference axis 


aerodynamic center defined by [ 0.45 - 


change in aerodynamic center due to flexibility, - X ac ^ ex j 
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z mcl 


coordinate of mean-camber line from wing reference plane, 


z upper + z lower 
2 


Z(x,y,t) force distribution in Z-direction (eq. (Ale)) 


a 



a. 



angle of attack measured from wing reference plane 
matrix defining angle of attack of each wing panel 
reference angle of attack 

matrix defining initial angle of attack at each wing panel due to camber, 
twist, or control deflection for wind off 


{ A “flex} 


matrix defining change in angle of attack at each panel due to aerodynamic 
loading 


/3 angle of sideslip 

6 deflection per unit load 


6 a wing trailing-edge control displacement in roll, positive for right aileron 

trailing edge down 

6 e wing trailing-edge control displacement in pitch, positive for trailing edge 

down 


6x ratio of modal deflection at station x to deflection at the station for which 

generalized mass is desired 


e 


v 

p 


pitch angle 

structural flexibility influence -coefficient matrix 

element of structural matrix [©) which defines the angle that panel j 
deflects due to a unit load on panel i 

nondimens ional spanwise coordinate, 

air density 
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frequency of ith structural mode 




Subscripts: 


A 

full-size airplane 

flex 

flexible model 

m 

model 

rigid 

rigid model 

sy 

symmetric 

as 

antisymmetric 

Matrix notation: 

1 1 

1 1 

square matrix 

LJ 

row matrix 

O 

column matrix 

□ 

square diagonal matrix 


Dots above symbols indicate derivatives with respect to time. 

DESCRIPTION OF MODELS 
Rigid Model 

The shape of the rigid model is dictated by performance and defined as the shape 
of the airplane "frozen" at the lg cruise condition (M = 2.7). This shape was predeter- 
mined by preliminary wind-tunnel tests at the NASA Langley Research Center and was 
used as a basis for all wind-tunnel tests during this investigation. A drawing of the com- 
plete model configuration is shown in figure 1. Detailed geometric properties of the 
model are presented in table I. Both models were built to a geometric scale factor of 
0.015. The solid steel (rigid) model incorporates a slender cambered body with a 74° 
swept -wing planform, four simulated flow-through engine nacelles, and two vertical tails 
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mounted outboard on the wing. The swept -wing planform has a subsonic leading edge at 
the cruise Mach number except in the region of the tip where the leading-edge sweep is 
reduced to 65°. 


Flexible Model 

The 0.015-size flexibly scaled model was designed from the shape of the rigid 
model at cruise condition (m = 2.7, Cl = 0.09, = 25.855 kN/m 2 ). This model was 

elastically scaled with a rigid forward fuselage and elastic wings and aft fuselage. The 
condition for model to aircraft similarity, according to reference 1 , is 



where (6m/® a) anc * (Shi/Sa) are the ratios of model to airplane deflection per unit load 
and reference wing area, respectively. A dynamic -pressure ratio (q m /q^J °f 1-20 
was arbitrarily selected. In order to simulate aerodynamic characteristics of the flex- 
ible airplane, it is necessary only to scale static load deflections; therefore, model weight 
was not scaled. A table of model scale factors is presented in table n. 

The flexible model was built so that at the cruise point the shape of the rigid and 
flexible models is the same. In order to achieve this condition, the flexible model is 
constructed to a model jig shape that is defined as the shape of the model when aerody- 
namic loads at the cruise point are removed. (The model jig shape will differ from that 
of the airplane jig shape since neither airplane mass nor mass distribution is simulated.) 
Figure 2 shows typical comparisons of the mean-camber line between cruise shape and 
the corresponding jig shape at various span stations. 

The flexible model was designed and constructed by The Boeing Company using a 
structural layout that would be quite similar to that of a full-size airplane. Structural 
ribs were fabricated from a balsa wood and fiber -glass sandwich with a thin aluminum 
cap. Plastic foam was used between structural spars to provide the proper wing contour. 
Finally, fiber -glass skins were bonded to the upper and lower surfaces to provide the 
properly scaled stiffnesses. Figure 3 shows a photograph of the flexible model during 
construction. A photograph of the complete flexible model prior to wind-tunnel testing 
is shown in figure 4. 


Flexible -Model Properties 

Because of the different structural representations required in the two analytical 
approaches, it was necessary to measure both a set of modal properties (mode shapes, 
frequencies, and generalized masses) and structural influence coefficients. Both sets of 
measurements were made for a sting -mounted model. 
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Modal Properties 


As input to the modal analysis, it is necessary to determine a set of generalized 
masses, mode shapes, and natural frequencies of model vibration. For this model, the 
first five symmetric structural modes, generalized masses, and natural frequencies were 
all that could be measured. 

Mode shapes and natural frequencies .- Vibration properties of the model were 
determined using the mode-shape measuring apparatus shown in figure 5. Two small 
electromagnetic shakers were positioned under the model to provide excitation for each 
of the first five symmetric modes. (Phasing between wing tips was checked for each 
mode to insure only symmetric response.) At each model frequency the displacement 
at 56 model control points was measured using a variable reluctance pickup. The loca- 
tions of these points on the wing and fuselage are shown in figure 6. A reference pickup, 
located near the wing tip, provided a continuous record of model displacement which is 
required to phase and correct raw data. Measured frequencies and nondimensional modal 
displacements (normalized to 1.0 at point 56) at the control points are presented in 
table HI. 

Generalized masses .- In order to characterize the mass properties of the model, 
it is not necessary to know the mass distribution explicitly because in the analysis the 
mass always occurs in the form of an integral value over the structure; that is, 

Mi = JJ m(x,y)hi 2 (x,y)dx dy 

This integral can be evaluated for each mode by using either the known mass distribution 
and mode shape or by experimental methods. Since the mass distribution was not deter- 
mined during construction, the direct experimental approach was used. 

The experimental technique employed is described in reference 4. A brief descrip- 
tion of this technique and its application follows: 

The generalized stiffness of the ith vibration mode is defined as 

Ki = MiCO.2 

where is the natural frequency of the ith structural mode. If it is assumed that by 
adding a small known mass Am, at station x, the generalized stiffness does not change, 
then 

Kj = (m 4 + AmSx 2 )^ 2 
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where 


6x ratio of modal deflection at station x to deflection at the station for which 

generalized mass is desired 

oJj natural frequency of ith mode with added mass Am 


Therefore, 

(Mi + Am 6x^) = 

or 


Mi- 


Am6x^ 



and 8x, a>j, and Wj are measured. The general- 
by the preceding equation. In practice, it is con- 

of — -1 for various values of Am in order 

\ w i/ 

to insure that the mode has not been altered by the additional mass. If the mode has been 
altered, this plot will be nonlinear. The slope of this curve, evaluated near zero for the 
nonlinear case, is the generalized mass. 


For each mode, Am is known 
ized mass Mj can now be evaluated 

venient to plot Amfix^ as a function 


The application of this technique to the flexible model is presented in figure 7. Also 
presented are the approximate locations of the symmetric node lines for each mode. For 
each mode a minimum of three data points were obtained. All generalized masses are 
referenced to model control point 56. 


Structural-Influence Coefficients 

The required structural information for the influence -coefficient analysis is the 
matrix jjsT| (eq. (B3)). This matrix defines the slope at any point on the structure due 
to a unit load at every other point. The flexibility influence coefficients for the flexible 
model were experimentally determined by The Boeing Company. 

The data required to define the flexibility influence coefficients were measured using 
the test setup shown in figure 8. The model was rigidly attached to a test fixture. Fifty - 
one electrical transducers were positioned over the right wing at preselected control 
points to monitor displacements. Figure 9 shows the paneling scheme and location of 
each control point. (Control points are at panel centroids.) The deflections at stations 1 
through 10 were not measured since the model was essentially rigid in this region. A 
pneumatic cylinder was used to apply point loads. 
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In order to evaluate the flexibility matrix, it is necessary to make two measure- 
ments: (1) the deflection on the right wing due to a load on the right wing, and (2) the 
deflection on the right wing due to a load on the left wing. These measurements were 
made by applying a known load at each of the 51 stations and monitoring the resulting 
deformation at all stations due to this load. If we define these two measured quantities as 


X 


matrix defining deflection on the right wing due to a unit load on the right wing 



matrix defining deflection on the right wing due to a unit load on the left wing 


then the symmetric deflection matrix can be expressed as 



and the antisymmetric deflection matrix as 



The matrix jjf] is now differentiated, by using a finite difference technique, to provide 
the required result [©]; that is, the slope at a point due to a unit load at any point on the 
structure. Measured values of deflection per unit load used to define matrices (aJ 
and (if) are given in table IV. 


WIND-TUNNEL METHODS AND TEST CONDITIONS 


The wind-tunnel studies were conducted in the Langley 8 -foot transonic pressure 
tunnel and in the Langley Unitary Plan wind tunnel in order to obtain data at subsonic and 
supersonic speeds. The stagnation dewpoint temperature was maintained sufficiently low 
to avoid any significant condensation effects in either tunnel. 

In order to insure a turbulent boundary -layer condition, all tests of both models 
were conducted with boundary -layer transition strips. These strips were composed of a 
band of No. 60 carborundum grit located 3.05 cm aft of the forebody apex and a band of 
No. 80 grit located 1.52 cm streamwise from the leading edge of all external surfaces 
and on the inside surface of the engine nacelles. 

Aerodynamic forces and moments about the 0.45c reference center were measured 
using a six-component strain-gage balance mounted within the models. Angle of attack 
was corrected for deflection of the sting and balance under aerodynamic load and for 
tunnel-flow angularity. Angle of attack was varied from about -10° to 5°, and sideslip 
angle was varied from -6° to 6°. 
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Wind-tunnel studies at M = 0.6, 0.9, and 1.2 over a range of dynamic pressures 
from 11.970 to 38.304 kN/m 2 were conducted in the Langley 8-foot transonic pressure 
tunnel. This tunnel is a variable -pressure, single -return facility having a slotted test 
section and is capable of producing velocities in the test section up to M = 1.3 without 
appreciable effects of choking and blockage. The nominal test conditions for this inves- 
tigation are: 


Mach 

number 

Stagnation 

temperature, 

°C 

Dynamic 

pressure, 

kN/m 2 

Reynolds number 
per meter 

Rigid model 

0.60 

48.9 

12.593 

6.562 X 10 6 

.90 

48.9 

17.093 

6.562 

1.20 

48.9 

19.966 

6.562 

Flexible model 

0.60 

48.9 

11.970 

6.201 X 10 6 

.60 



23.940 

12.434 

.60 



32.654 

16.962 

.90 



11.970 

4.593 

.90 



23.940 

9.186 

.90 



29.638 

11.385 

1.20 



11.970 

3.937 

1.20 



23.940 

7.874 

1.20 



38.304 

12.598 


Wind-tunnel studies at M = 2.3 and M = 2.7 over a range of dynamic pressures 
from 11.970 to 47.880 kN/m^ were conducted in the Langley Unitary P lan wind tunnel. 
The nominal test conditions for this investigation are: 


Mach 

number 

Stagnation 

temperature, 

°C 

Dynamic 

pressure, 

kN/m 2 

Reynolds number 
per meter 

Rigid model 

2.3 

65.6 

21.738 

6.562 X 10 6 

2.7 

65.6 

19.822 

6.562 

Flexible model 

2.3 

65.6 

11.970 

3.642 X 10 6 

2.3 



23.940 

7.251 

2.7 



11.970 

3.970 

2.7 



23.940 

7.940 

2.7 



31.122 

10.335 

2.7 



47.880 

15.879 
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AEROELASTIC ANALYSIS 


Analytical calculations for the rigid model as a function of Mach number and for 
the flexible model as a function of Mach number and dynamic pressure have been made 
employing the modal method (appendix A) and structural influence coefficients (appen- 
dix B). The calculations were made for a sting -mounted, "massless" model. The term 
massless is used to describe the case in which only aerodynamic loads are considered 
acting upon the structure; that is, inertia loads such as weight times load factor are 
assumed to be zero. 


Modal Method 

A digital computer program, using the approach discussed in appendix A, was 
developed to predict aeroelastic effects on longitudinal stability and control characteris- 
tics. Calculations were made for a sting-mounted model using the first five measured 
symmetric modes and generalized masses (see fig. 7 and table in). 

The basic free -flight equations presented in appendix A represent the sting -mounted 
model, if it is assumed that integrals of the form 

JJ m(x,y)hi(x,y)dx dy 

are negligibly small. In the free -flight case these integrals are exactly zero, except for 
the plunge mode where the integral simplifies to the total mass of the airplane, due to the 
orthogonality of the structural and rigid-body modes. 

Subsonic aerodynamics . The subsonic generalized aerodynamic forces 
JJ Ap i (x,y)hj(x,y)dx dy 

were formulated through the use of kernel -function aerodynamics. The pressure dis- 
tribution Apj(x,y) is obtained by a numerical method similar to that presented in refer- 
ence 5. The pressure distribution is obtained by numerically solving a linear integral 
equation which relates the pressure distribution to the downwash. 

The solution of the kernel equation for Ap^(x,y) requires that the downwash be 
specified at known control points. Thirty -six control points were used in these calcula- 
tions. These points were located at 12.1, 33.6, 56.7, 76.5, 88.5, and 96.8 percent of the 
panel span ( 77 ) and at 5.7, 21.6, 44.0, 67.7, 87.4, and 98.5 percent of the local streamwise 
chord (x/c). Only the wing, extended to the fuselage center line, is represented 
aerodynamically. 
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Substituting the appropriate expressions for the downwash due to rigid loading and 
that due to flexibility, the generalized aerodynamic forces required to solve the equa- 
tions of motion (appendix A) are determined. 

Supersonic aerodynamics .- Supersonic generalized aerodynamic forces were evalu- 
ated through the use of a supersonic Mach-box procedure (ref. 6). As in the kernel pro- 
gram the pressure distribution is obtained by specifying the downwash at known control 
points. The 56 control points at which the flexible modes were specified (fig. 6) were 
also used as the downwash control points for this calculation. Only the wing, extended 
to the fuselage center line, was represented aerodynamically. 

Once the pressure distribution due to rigid loading and each flexible mode are 
known, the generalized aerodynamic forces are calculated. The equations of motion 
(appendix A) are now solved as before. 

Influence -Coefficient Method 

A digital computer program developed by The Boeing Company, using the approach 
discussed in appendix B, was used to predict aeroelastic effects on longitudinal stability 
and control and lateral control characteristics. Calculations were made for a sting - 
mounted model using measured flexibility influence coefficients (see table IV). 

The free -flight equations presented in appendix B represent the massless, con- 
strained model if the weight matrix (w) is assumed to be 0. Equation (B4), therefore, 
becomes 

{F} = q[ S jA]0 

When the panel weights are set equal to zero, the aeroelastic characteristics become a 
function of planform and structure at a given test condition. Comparisons of aerody- 
namic characteristics with and without mass effects are presented in table V. 

The aerodynamic matrix [a] was formulated through the use of constant -pressure 
panels that can be used at subsonic and supersonic speeds (ref. 7) and which were coinci- 
dent with the structural panels (fig. 9). The aerodynamic representation included the 
forebody, wing, and cylindrical afterbody. 

Once the area matrix , aerodynamic matrix [a] , and flexibility influence - 
coefficient matrix [©] are evaluated, the equations in appendix B are solved for the 
flexible aerodynamic characteristics as a function of rigid loading and test condition. 

RESULTS 

This section compares longitudinal aerodynamic characteristics obtained with the 
rigid model and the flexible model and from theoretical calculations. Also included is a 
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presentation of lateral experimental results obtained for the flexible model and a com- 
parison of lateral-control effectiveness with influence -coefficient calculations. The 
rigid model is assumed perfectly rigid, and no corrections to the data due to flexibility 
are included. 


Longitudinal Aerodynamic Characteristics 

A comparison of analyses and wind-tunnel data in predicting the variation of lift- 
curve slope CL a and aerodynamic center X ac , in percent of the reference chord, as 
a function of Mach number is presented in figure 10. Measured results are given for 
both the rigid model (assumed independent of dynamic pressure) and the flexible model 
at a dynamic pressure 23.94 kN/m2. The effect of dynamic pressure at all Mach num- 
bers is to reduce CL a significantly and to shift the aerodynamic center in a destabiliz- 
ing forward direction. Reductions in Cj_, a of about 37 percent and shifts in the aero- 
dynamic center as large as 18 percent of c are evident at M = 1.2. These phenomena 
are associated with structural deflections outboard on the wing which tend to reduce the 
local angle of attack and, subsequently, to reduce the contribution of the outboard portion 
of the wing to the total lift. This reduction in lift results in a lower lift -curve slope and 
a forward shift in the aerodynamic center. In general, correlation between analyses and 
experiment is reasonable. 

The effect of dynamic pressure on lift-curve slope is predicted somewhat better 
than the effect of dynamic pressure on aerodynamic-center location. Figure 11 com- 
pares analyses and wind-tunnel data in predicting the variation of the ratio of flexible 
C L a to rigid CLa and aerodynamic -center movement between rigid and flexible 
models as a function of dynamic pressure. Data are presented for Mach numbers of 0.6, 
0.9, 1.2, 2.3, and 2.7 over a range of dynamic pressure. The shift in aerodynamic center, 

AX ac = (x ac>r igid - X ac i ex ) 

is a measure of the movement of the aerodynamic center from its location on the rigid 
model as indicated in figure 10. In general, the agreement between analyses and wind- 
tunnel data in predicting the variation in CL a is good. For most cases, the agree- 
ment between modal and influence -coefficient techniques is comparable; however, the 
aerodynamic -center movement with dynamic pressure is predicted better with the modal 
technique at M = 0.6, 0.9, and 1.2. Both techniques show good agreement with wind- 
tunnel data at M = 2.3 and M = 2.7. It should be noted that at these higher Mach num- 
bers the shift in aerodynamic center is predicted quite well even though neither analysis 
predicted the absolute position of the aerodynamic center (see fig. 10) . 

Figure 12 presents a comparison of analyses and experiment in predicting lift and 
pitching -moment coefficients at a = 0° as a function of dynamic pressure. The cal- 
culated data are generated by putting into the program the camber and twist distribution 
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of the wing mean-camber line (fig. 2). Only flexible -model results are presented since 
the slope of the mean-camber line for the rigid model applies only to the design cruise 
point. It should be noted that calculated data at a dynamic pressure of zero correspond 
to a rigid configuration with the flexible -model jig-shape camber and twist. 

At Mach numbers of 0.6, 0.9, and 1.2, the influence -coefficient analysis appears to 
predict the pitching -moment coefficients somewhat better than the modal approach. At 
Mach numbers of 0.6 and 0.9, both methods are directly comparable for lift coefficient. 

At M = 1.2, the influence -coefficient method predicts both lift coefficient and pitching- 
moment coefficient somewhat better. The results at M = 2.3 and 2.7 indicate that both 
methods predicted the lift coefficient quite well, but for these cases the modal method 
appears to predict the pitching -moment coefficient somewhat better. 

The reason one method appears to give better results at one Mach number than at 
others is probably the aerodynamic representations utilized in both methods. The modal 
method uses kernel function subsonically and Mach-box aerodynamics supersonically, 
whereas the influence -coefficient technique uses constant -pressure panels throughout 
and includes fuselage effects. In all cases the trend is accurately predicted, and in most 
cases the increment between measured data is estimated quite well using either technique. 

A comparison between wind-tunnel data and analyses in predicting control-surface 
effectiveness in pitch is presented in figure 13. At M = 0.6 and 0.9 only influence - 
coefficient analysis data are presented since the modal technique could only be applied 
to the supersonic case. (The kernel-function program used does not provide for inclu- 
sion of control-surface aerodynamics.) Data at Mach numbers of 1.2, 2.3, and 2.7 com- 
pare both influence -coefficient and modal analyses with experimental results. Wind- 
tunnel data for the rigid model were not available at these Mach numbers for the control 
surface under investigation. 

Correlation between experiment and analyses is quite poor for predicting the 
absolute level. The increment in control effectiveness, however, is predicted somewhat 
better. It is believed that part of this problem is due to the inadequacy of both aerody- 
namic theories (Mach-box and constant -pressure panels) to properly predict the rigid 
loading because of the proximity of the control surface to both the vertical tail and the 
outboard engine. The problem was further aggravated by the flexibility of the control- 
surface attachment points. Since both analytical methods at M = 1.2, 2.3, and 2.7 pre- 
dict about the same rigid loading, the differences in calculations are attributed to the 
limited number of modes used in the modal calculation. 

Lateral Aerodynamic Characteristics 

The variation of measured lateral stability derivatives C^, C n ^, and Cy^ for 
the flexible model as a function of angle of attack at different dynamic pressures is 
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presented in figure 14. These derivatives were estimated around /3 = 0° by taking, 
for example, 



( C ^/3=2° ~ 
2 


At present the theoretical calculations have not been extended to predict these lateral- 
stability characteristics. The general analytical approach outlined previously can be 
applied to the lateral case by using the appropriate antisymmetric structural and aero- 
dynamic representations. 

For the Mach number and dynamic-pressure range investigated, the flexible model 
exhibited positive effective dihedral (-Q^) except at the larger negative angles of attack, 
a positive value of C n ^ except at the larger negative angles of attack at M = 2.3 and 2.7, 
and a negative side-force derivative Cy^. For this configuration, dynamic pressure 
appears to have a less significant effect on the lateral -directional characteristics than 
on the longitudinal characteristics. The primary effect is a reduction in C at all 
Mach numbers, with increasing dynamic pressure. 

The effect of dynamic pressure on control -surface effectiveness in roll is pre- 
sented in figure 15. Calculations using the influence -coefficient technique were made to 
predict this effect. At all Mach numbers the effect of an increase in dynamic pressure 
is to reduce appreciably the aileron effectiveness in roll. At M = 1.2 and q = 33 kN/m2 
the model exhibits zero control effectiveness. Once again, the correlation between anal- 
ysis and experiment in predicting control-surface effectiveness is quite poor. As stated 
previously, this can be attributed to poor prediction of the rigid loading. It is interest- 
ing to note that zero effectiveness in pitch (M * 1.2, fig. 13) and aileron reversal in roll 
(M = 1.2, fig. 15) are analytically predicted reasonably close to the experimental results 
using the influence -coefficient analysis. 


CONCLUDING REMARKS 


Both wind-tunnel studies and analytical calculation of 0.015-size rigid and flexibly 
scaled models of a proposed supersonic transport configuration have been conducted at 
subsonic and supersonic speeds to measure the effect of flexibility on the aerodynamic 
characteristics. Analytical calculations using both a modal approach and structural 
influence -coefficient approach are presented. 

Examination of the analytical and experimental data indicates that: 

1. The analyses did predict reasonably well the effect of flexibility on the basic 
longitudinal characteristics; however, the analyses are shown to be poor in predicting 
control-surface derivatives in both pitch and roll. 
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2. Both the modal approach and the structural influence -coefficient approach yield 
generally comparable results. 

Even though the primary purpose of this investigation was to evaluate analytical 
techniques, it has been shown experimentally that: 

1. The longitudinal aerodynamic characteristics of the flexible model were affected 
strongly by flexibility effects which included large reductions in lift -curve slope, desta- 
bilizing shifts in the aerodynamic center, and large reductions in control effectiveness. 

2. The lateral aerodynamic characteristics of the flexible model were not greatly 
affected by flexibility except for control effectiveness in roll. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 10, 1972. 
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APPENDIX A 


MODAL METHOD FOR DETERMINING LONGITUDINAL 
STABILITY AND CONTROL CHARACTERISTICS 

A detailed analysis of the procedure for calculating aeroelastic effects on longi- 
tudinal stability and control using the modal approach is presented. 

The deformation of the airplane as a function of rigid loading is determined by 
solving Lagrange’s equations of motion. From these calculated deformations the effec- 
tive wing loading is obtained which, in turn, is used to calculate the flexible stability and 
control characteristics. 

A digital computer program, using the approach discussed in reference 8, was 
developed to predict aeroelastic effects on longitudinal stability and control. 

Equations of motion .- Under the following assumptions 

(1) The forward speed V of the airplane is constant 

(2) Small angle -of -attack variation ~ Cj_J 

(3) Small structural deformation 

(4) Structural motion is slow enough that structural acceleration and 

structural velocity are negligibly small 

(5) Orthogonal structural modes 

the linearized free -flight longitudinal equations of motion appear in a general form (in 


the body -axis system) as follows: 

M i q i (t) + w i 2 M i q i (t) = Qj(t) (Al) 

where 

M i = m(x,y)h i 2 (x,y)dx dy (Ala) 

Qi(t) = Z(x,y,t)h i (x,y)dx dy (Alb) 

Z(x,y,t) = Ap(x,y,t) - ngm(x,y) (Ale) 

m 

Ap(x,y ,t) = ^ APj(x,y,t)q.(t) (Aid) 

3=1 
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APPENDIX A — Continued 


Generalized forces and moments .- In order to determine the generalized forces 
and moments Qj(t), it is necessary to relate the downwash w(x,y,t) to the displace- 
ment H(x,y,t) of the system. It is assumed that the displacement H(x,y,t) can be 
represented by a superposition of the normal modes of vibration so that 

H(x,y,t) = hjfoyjq^t) + h 2 (x,y)q 2 (t) + . . . + h m (x,y)q m (t) 

where q^ specifies the magnitude of the displacement in the ith mode, and hi(x,y) 
gives the shape of the mode. Letting hj(x,y) be rigid-body plunge and h 2 (x,y) be 
rigid-body pitch, then 

hi(x,y) = 1 h 2 (x,y) = x 

qj(t) = z q 2 (t) = 6 

The displacement can now be written as 

m 

H(x,y,t) = z + x6 + ^ hi(x,y)q.(t) 
i=3 

The downwash w(x,y,t) associated with the displacement H(x,y,t) is given by 


w(x,y,t) = (V ■^ + ^) H (x,y,t) 


Downwash due to angle of attack .- Assuming the 
small and using assumption (4), q^ = 0 where i = 3, 
may be expressed as 

m 8 hi(x,y) 


w(x,y) = V 


I 

i=3 


9x 


q A + a 


downwash associated with 6 is 
the downwash due to angle of attack 


(A2) 


where 



(A2a) 


The downwash is composed of two parts; namely, 

V ah i(x,y) 

Downwash due to flexibility: V > q, 

6 * 1 

Downwash due to rigid loading: V a 

Once the downwash distribution is known over the wing (eq. (A2)), the pressure dis- 
tribution Ap i for each flexible and rigid mode can be obtained from available aerody- 
namic theories. In the modal framework, pressure distributions are determined in the 
following manner: 
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APPENDIX A — Continued 



Determination of structural deformation .- Using free -free normal modes (both 
rigid-body and structural) and downwash due to angle of attack, the longitudinal equa- 
tions of motion (Al) can be written as 

MZ = q 3 JJ Ap 3 (x,y)dx dy + . . . + q m 

+ a Ap a (x,y)dx dy - Mgn (A3) 

15 - <13 II x Ap 3 (x,y)dx dy + . . . + q m x Ap m (x,y)dx dy 

+ OL JJ X Ap a (x,y)dx dy (A4) 


yy APmfrtf)* 1 * dy 


a '3 2M 3 <1 3 = q 3 II h 3< x >y) ^P 3 (x,y)dx dy + . . . + q m Jj* hg(x,y) Ap m (x,y)dx dy + a hg(x,y) Ap a (x,y)dx dy 

l 

"m 2M m q m = q 3 II Ap 3 (x,y)dx dy + . . . + q m H A Pm(*>y)dx dy + a h^x.y) Ap a (x,y)dx dy 

where 

a>i = cog = 0 

M = y^ m(x,y)hj 2 (x,y)dx dy = §§ m(x,y)dx dy 

I = yy m(x,y)h 2 2 (x,y)dx dy = §§ x 2 m(x,y)dx dy 
and integrals of the form 

ng yy m(x,y)hi(x,y)dx dy = 0 (i = 2, 3, . . m) 

due to orthogonality of the free -free modes. 


(A5) 
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APPENDIX A - Concluded 


Once the generalized masses Mj, the natural structural frequencies and the 

aerodynamic forces per unit generalized displacement Cj* h^(x,y) Apj(x,y)dx dy are 

evaluated, equations (A5) are solved simultaneously for q i? the generalized displace- 
ment q^ per unit angle of attack a. 

Final form of flexible characteristics .- Defining normal force and pitching moment 
per unit generalized displacement as 

N i = JJ Ap i (x,y)dx dy ^ x Ap i (x,y)dx dy 

the plunge equation (A3) can be expressed in terms of q^ as 


MZ = a(q 3 N3 + q 4 N 4 + 


• • • + q m N m + 


N a ) - Mgn 


(A6) 


The aerodynamic terms in equation (A6) can be equated to the normal -force coefficient as 

|pV 2 SC N = «(q 3 N 3 + q 4 N 4 + . . . + q m N m + N a ) 

Differentiating with respect to a, the flexible lift -curve slope is obtained as a function 
of the modal properties q i and N| as 

C k a;, flex * C Nq! } f lex = I~2^(^3 N 3 + ^4 N 4 + • • • + V N m) + C N a>rigid < A7 ) 

Using the pitching -moment equation (A4), the static stability derivative C mQ! ^ can 
be expressed in terms of and Mj as 


Cm o', flex “ T~^( q 3 M 3 + ^4 M 4 + • • • + q m M m) + c m Q , jrigid 


(A8) 


In a similar manner, the aerodynamic coefficients, such as C L 0 , C m Q , , 

9 9 

C m6e , C L ^, and C m q can be evaluated by using the appropriate loading conditions, 
that is, by specifying the proper downwash distribution associated with each rigid loading. 
In order to evaluate c L,Ofi ex 311(1 C mj0 ^ ex , for example, replace integrals of the form 

a JJ hj( x ,y) Ap a (x,y)dx dy 
in equations (A3) to (A5) with integrals of the form 

jl M x .y) A P( x >y) c araber and twist * dy 


and solve as before. 
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APPENDIX B 


INFLUENCE -COEFFICIENT METHOD FOR DETERMINING LONGITUDINAL 
STABILITY AND CONTROL CHARACTERISTICS 


A detailed analysis of the procedure developed by The Boeing Company for cal- 
culating the aeroelastic effects on longitudinal stability and control characteristics using 
the influence -coefficient approach is presented. 

The airplane is first divided into a number of panels (see fig. 9). For each panel 
an aerodynamic matrix element A ij is developed which relates the pressure at panel i 
due to a unit rotation at panel j. A structural flexibility matrix element 0jj is obtained 
which relates a change in angle of attack at panel j due to a unit load at panel i. With 
these matrices, which define the aerodynamics and flexibility of the airplane, the follow- 
ing set of matrix equations can be written: 


( L ) = «[ s ] [A] 

(Bl) 

{»} = “r© + {“q> + (Aaflex) 

(B2) 

{ a “nex}= 

(B3) 

{F} = q[s][A]{ a }-n{w} 

(B4) 


In order to solve for lift-curve slope, the preceding equations can be solved in the fol- 
lowing manner: 

= «r® + {“q) + [ 0 ]( F > 

Substitute for {f} , equation (B4); then 

{«} = “r® + <S) + [ 0 ](q[ s ][ A ]{ a } - n (w» 

Solving for ^oTj> and substituting this result into equation (Bl) results in the total lift 

L = LlJ {L> - «LlJH [A] {(H - «HEI[a]) '‘(air© ♦ <j>o) - n[e] (w»j 


Since the load factor 
expressed as 


n 


is defined as 


Lift 

Weight 


and 


T i -ft- 

C L = the previous equation can be 


C L - ~ ((H - 4[®][ s ][A])‘ l (“r{0 + {<> 
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APPENDIX B - Concluded 


Differentiating this expression with respect to a T and simplifying results in the final 
form of the lift -curve slope; that is, 

Cl Marais (B5) 

s(i ♦ |L1JH[A](H - q [eXs][A] )- 1 Qe] <w>) 

Since the lift Lj on each panel is known, the pitching -moment derivatives can be readily 
calculated from panel geometry; that is, 

In a similar manner, expressions can be derived for other aerodynamic characteristics 
such as Cl j0 , C m>0 , Oj^, Cm§ e > ^Lq> an< * ^iHq* 

The principal advantage of the direct influence -coefficient approach is the capability 
for separation of inertia loading effects from aerodynamic loading effects. A comparison 
of these characteristics is presented in table V. The aerodynamic characteristics pre- 
sented are separated into three categories: aeroelastic effects due to both flexibility and 
inertia, aeroelastic effects due only to flexibility (referred to as massless derivatives in 
ref. 1), and rigid -airplane characteristics. 
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TABLE I 


GEOMETRIC PROPERTIES OF MODELS 


Wing 

Aspect ratio 1.630 

Span, meters 0.580 

Area, meters^ 0.207 

Root chord at fuselage center line, meters . . 0.815 

Tip chord, meters 0.034 

Reference chord, meters 0.487 

Fuselage 

Length, meters 1.167 

Diameter at base, meters 0.045 

Vertical tails 

Area, meters^ 0.005 

Thickness -chord ratio 0.029 


TABLE n 

MODEL SCALE FACTORS 


Quantity 


Formula 

Factor 

Length 

L m/ L A 

Selected 

0.015 

Dynamic pressure 

V/^A 

©G0* seiecied 

1.2 

Mach number 

M m /M A 


1.0 



/n \ /-r 


Air loads 

Fm/FA 

Wte) 

2.7 X 10 -4 

Deflection per unit load 

W«A 

gfe) 

55.6 

Scaling rule 

. /Vi/V 

' vWW, 

,/s^x 1 / 2 

y =i ° 














TABLE m 


[ONAL MODE -SHAPE DATA 


7 ] = 0.567 
c = 0.2933 m 

77 = 0.765 
c = 0.1430 m 

7 ] = 0.885 
c = 0.0889 m 

V = 0.968 
c = 0.0498 ra 

x/c 


x/c 

hi 

x/c 

hi 

x/c 

hi 

a) * 27.19 hertz, 

isl 







0.132 

0.055 

0.401 

0.054 

0.580 

0.061 

0.854 

.171 

.155 

.183 

.447 

.180 

.622 

.193 

.873 

.304 

.201 

.312 

.484 

.303 

.654 

.315 

.896 

.437 

.242 

.441 

.519 

.429 

.682 

.447 

.916 

.569 

.300 

.571 

.562 

.554 

.718 

.563 

.922 

.702 

.360 

.700 

.606 

.680 

.735 

.695 

.927 

.825 

.436 

.830 

.651 

.811 

.777 

.822 

.939 

.967 

.518 

.965 

.722 

.940 

.802 

.954 

1.000 

a ) - 42.87 hertz, 

i = 2 






0.038 


0.055 

0.216 

0.054 

0.548 

0.061 

0.780 

.171 


.183 

.272 

.180 

.563 

.193 

.799 

.304 


.312 

.329 

.303 

.603 

.315 

.840 

.437 


.441 

.393 

.429 

.653 

.447 

.860 

.569 


.571 

.467 

.554 

.691 

.563 

.884 

.702 


.700 

.530 

.680 

.741 

.695 

.917 

.825 


.830 

.603 

.811 

.805 

.822 

.958 

.967 

.417 

.965 

.689 

.940 

.855 

.954 

1.000 

Co = 68.21 hertz, 

i = 3 






0.038 

-0.001 

0.055 

0.090 

0.054 

0.394 

0.061 

0.693 

.171 

-.002 

.183 

.134 

.180 

.441 

.193 

.732 

.304 

-.003 

.312 

.177 

.303 

.486 

.315 

.783 

.437 

-.002 

.441 

.226 

.429 

.531 

.447 

.803 

.569 

-.001 

.571 

.285 

.554 

.582 

.563 

.846 

.702 

.005 

.700 

.336 

.680 

.642 

.695 

.891 

.825 

.044 

.830 

.401 

.811 

.702 

.822 

.933 

.967 

.116 

.965 

.462 

.940 

.745 

.954 

1.000 

a) = 103.09 hertz, 

i = 4 






0.038 


0.055 

0.046 

0.054 

0.280 

0.061 

0.671 

.171 


.183 

.056 

.180 

.317 

.193 

.683 

.304 


.312 

.071 

.303 

.351 

.315 

.722 

.437 


.441 

.086 

.429 

.394 

.447 

.753 

.569 


.571 

.109 

.554 

.443 

.563 

.809 

.702 


.700 

.130 

.680 

.483 

.695 

.854 

.825 


.830 

.158 

.811 

.551 

.822 

.932 

.967 


.965 

.188 

.940 

.618 

.954 

1.000 

w = 148.65 hertz 

, i * 5 






0.038 

-0.160 

0.055 

-0.809 

0.054 

-0.372 

0.061 

0.164 

.171 

-.308 

.183 

-.646 

.180 

-.263 

.193 

.280 

.304 

-.477 

.312 

-.480 

.303 

-.126 

.315 

.389 

.437 

-.540 

.441 

-.298 

.429 

-.001 

.447 

.483 

.569 

-.492 

.571 

-.112 

.554 

.077 

.563 

.610 

.702 

-.284 

.700 

-.005 

.680 

.257 

.695 

.725 

.825 

-.032 

.830 

.148 1 

.811 

.461 

.833 

.868 

.967 

.201 

.965 

.345 

.940 

.644 

.954 

1.000 
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DO 

00 


TABLE IV 

STRUCTURAL INFLUENCE -COEFFICIENT MEASUREMENTS 
(a) Matrix [a] 


Panel deflection per unit load, metera/newton, measured at position: 


point 

11 

12 

13 

14 

15 

16 

17 

IB 

19 

20 

21 

, 22 

23 

24 

25 

26 

27 

11 

.O37X10 -4 

• C46X10" 4 

.057X10 -4 

.009X10 ” 4 

.012X10" 4 

■ 017xl0“ 4 

.022XKT 4 

.02 9X10 " 4 

•039X10 -4 

.0 50X10 " 4 

•Q58X1Q -4 

•061X10 -4 

.015X10” 4 

.oiaxio -4 

.Q23X1Q-* 

*030X10 " 4 

.038X10 

12 

• 051 

• C < 9 

.107 

.008 

.010 

.020 

.027 

.037 

.058 

.084 

.103 

.116 

.018 

.021 

.027 

.030 

.054 

13 

• 071 

.118 

.205 

.010 

.015 

.022 

.030 

• 048 

.078 

.126 

.171 

.198 

.019 

• 023 

.029 1 

.047 

.071 

14 

.006 

.CC6 

.006 

• 014 

.006 

.006 

.007 

.007 

.007 

.007 

.007 

.007 

.000 

•00 8 

.009 

.009 

.014 

15 

.013 

• C 13 

.013 

.009 

.016 

.012 

.013 

• 014 

.014 

.012 

.015 

• CL5 

.018 

.016 

.017 

•017 

.017 

16 

• 018 

.019 

.021 

.010 

•013 

.025 

.021 

• 023 

.024 

.026 

• 026 

• 026 

.023 

• 029 

• 026 

.027 

.027 

17 

• 029 

.029 

.031 

.011 

.CIS 

.021 

.036 

.029 

.031 

.034 

.035 

.035 

.022 

.027 

.034 

.034 

.034 

18 

.037 

.040 

.045 

.012 

.016 

.023 

.029 

.038 

.041 

.046 

• 049 

.050 

.023 

.027 

.036 

.044 

.047 

19 

.045 

.0*5 

• 068 

.012 

.015 

.022 

.026 

.035 

.054 

.065 

.074 

.079 

.021 

•024 

.021 

.043 

.057 

20 

.059 

.067 

.119 

.011 

.016 

.025 

.031 

.044 

.070 

.083 

.130 

.146 

.023 

.026 

.034 

.049 

• 072 

21 

.071 

.114 

• 166 

.013 

.017 

.026 

.034 

.049 

.083 

.135 

.246 

.335 

.025 

•028 

.037 

.056 

.085 

22 

.063 

.111 

.173 

.011 

.014 

.021 

.027 

.041 

,076 

.135 

.312 

.678 

.020 

.022 

.029 

.045 

.075 

23 

.018 

.018 

.019 

.014 

.019 

.023 

.022 

.023 

.024 

.025 

.026 

.025 

.083 

.039 

• 032 

.031 

.031 

24 

• 023 

.023 

.02* 

.017 

.018 

.028 

• 026 

.027 

.028 

.030 

.031 

.030 

.040 

.066 

• 042 

.037 

.035 

25 

.033 

.033 

.034 

.017 

.020 

.027 

.035 

.038 

.039 

.041 

.041 

•040 

.022 

• 043 

• 065 

.054 

.050 

26 

.039 

• 042 

.049 

• 020 

.018 

.025 

.033 

.042 

.048 

.053 

.054 

.054 

.027 

• 029 

• 048 

.063 

• 062 

27 

• C49 

.c;e 

.069 

.016 

.019 

.026 

.032 

.045 

.062 

.074 

.082 

.082 

.028 

.033 

.044 

.061 

• 083 

28 

• 059 

• C 61 

.102 

.015 

.019 

.027 

.035 

.038 

.076 

.104 

.123 

.133 

.029 

.034 

.046 

.066 

.092 

29 

• 068 

* 1C3 

- i3e 

.015 

.019 

.029 

.037 

.053 

.086 

.133 

.180 

.215 

.029 

.033 

• 045 

.067 

.100 

30 

• 070 

.114 

.160 

.014 

.018 

.027 

.035 

• 051 

.090 

.148 

.230 

• 309 

.027 

.031 

.041 

•063 

.100 

31 

.065 

.113 

.167 

.010 

.015 

.021 

.029 

.043 

.085 

.152 

• 256 

.348 

.022 

• 025 

.022 

.054 

.092 

32 

.027 

• C 2 7 

.027 

.017 

.020 

.029 

.034 

.038 

.030 

.038 

.037 

.034 

.037 

.056 

.073 

.059 

.053 

33 

.035 

.035 

.035 

.019 

.021 

.029 

.037 

.044 

.046 

.048 

.047 

• C45 

.036 

.046 

.070 

.073 

.066 

14 

.046 

• C so 

.053 

.018 

.022 

.030 

.038 

.051 

.062 

.068 

.070 

.069 

.036 

.04 3 

.062 

.082 

.090 

35 

.057 

.069 

♦ 079 

.019 

.023 

.032 

.041 

.055 

.078 

.095 

• 096 

.107 

.037 

.043 

.059 

.083 

.110 

36 

.067 

• C 69 

• 109 

.020 

.025 

.034 

• 042 

.059 

.091 

.122 

.144 

.157 

.039 

.043 

.059 

.085 

.119 

37 

.074 

, *K7 

.139 

.020 

.025 

.035 

.044 

.063 

.100 

.146 

.188 

.219 

.039 

.043 

.059 

.087 

.127 

38 

.079 

.123 

.164 

.019 

.025 

.034 

.043 

.063 

.106 

.165 

.234 

.290 

.029 

.042 

.057 

.086 

.130 

39 

.082 

.121 

.182 

.018 

.023 

.033 | 

.042 

.062 

.110 1 

.179 

.269 

.344 

.036 

.039 

.053 

• 083 

• 131 

40 

.069 

.120 

.174 

.011 

.016 

.024 

.030 

.047 

.095 

.170 

.273 

.360 

.025 

.023 

.037 

.065 

.113 

41 

.059 

• C<9 

.077 

• 024 

.028 

.036 

.044 

• 061 

.085 

.103 

• 110 

• 110 

.046 

.050 

.071 

.100 

.131 

42 

.063 

• C 76 

.009 

.023 

.028 

.037 

.046 

.064 

.093 

.115 

.127 

.131 

.047 

.051 

.071 

.102 

• 138 

43 

• 0 65 

.Cf7 

.104 

.022 

.025 

.034 

.042 

• 060 

.095 

.128 

.149 

• 160 

.041 

.045 

• 063 

.096 

.137 

44 

• C67 

.097 

.120 

.021 

.025 

.033 

.040 

.059 

.099 

.138 

.176 

.200 

.02E 

.04 2 

• 060 

.093 

.139 

45 

• C70 

• 106 

• 138 

.022 

.023 

.031 

.039 

.059 

.102 

.154 

.204 

.24 2 

.026 

• 038 

• 055 

.089 

.138 

46 

.074 

.115 

.154 

.017 

.038 

.030 

.038 

.058 

.106 

.167 

.235 

.288 

.024 

.037 

• 053 

.087 

• 140 

47 

.076 

.126 

.172 

.018 

.022 

.030 

.037 

.058 

.111 

.182 

.266 

.337 

.032 

.034 

.051 

.086 

.142 

48 

.071 

.128 

.181 

.016 

.018 

.023 

.030 

.051 

.105 

.185 

.286 

.371 

.025 

•026 

.040 

.074 

.132 

49 

.073 

.134 

.192 

.015 

.018 

.025 

.030 

.050 

.108 

.194 

• 308 

.404 

.025 

•025 

.040 

.074 

<135 

50 

.069 

• 112 

.140 

.017 

.019 

.025 

.032 

.054 

.104 

.167 

.227 

• 271 

.026 

.026 

.048 

.089 

.150 

51 

.070 

• 115 

.151 

• 016 

.019 

.024 

.033 

.054 

.108 

.175 

.246 

.297 

*027 

.02 8 

.046 

.089 

• 151 

52 

.074 

• 122 

.162 

.014 

.017 

.025 

.033 

.057 

• 114 

.184 

.265 

.328 

.025 

.027 

.048 

.090 

.156 

53 

.077 

• 128 

.174 

.015 

.018 

.025 

.034 

.056 

.116 

.192 

.283 

.354 

.027 

.027 

.047 

.087 

.155 

54 

.078 

.132 

.179 

.014 

.017 

.024 

.033 

.057 

• 116 

.195 

.291 

.371 

.026 

.026 

.047 

.086 

.158 

55 

.079 

.140 

• 190 

.014 

.018 

.026 

.035 

.059 

.121 

.209 

.317 

.404 

.026 

.027 

.040 

.091 

.163 

56 

.079 

.144 

.202 

.018 

•021 

.026 

.032 

.056 

.123 

.214 

.334 

• 431 

.027 

.028 

.043 

.087 

.157 

57 

.085 

.1*0 

.208 

.015 

.018 

.025 

.034 

.059 

.130 

.224 

.340 

.452 

.027 

.025 

.047 

.C90 

.165 

58 

.082 

.142 

.191 

.013 

.018 

.026 

.035 

.062 

.130 

.218 

• 321 

.405 

.023 

.026 

.050 

.099 

.180 

59 

.064 

.148 

• 204 

.016 

.019 

.025 

.034 

.061 

.132 

.227 

.340 

.437 

.025 

.024 

.049 

.096 

.179 

60 

.oe6 

.155 

.213 

.013 

.018 

.026 

.033 

.059 

.138 

.234 

.356 

.459 

.023 

•026 

,047 

.097 

*180 

61 

.086 

.159 

.223 

.014 

.017 

.025 

• 033 

.061 

.137 

.242 

.376 

.489 

.023 

.024 

.047 

.097 

.182 
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TABLE IV 

STRUCTURAL INFLUENCE-COEFFICIENT MEASUREMENTS - Continued 
(a) Matrix [aT] - Concluded 


Load 

Panel deflection per unit load, meters/newton, measured at position: j 


point 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

.60 

"61 


11 

•CC2X10" 4 

.04 8X10 " 4 

*05 5X10 ” 4 

.057X1 CT 4 

• 050X10 -4 

• 046X10 -4 

.057X10" 4 

. 054X10 " 4 

• 066X10 -4 

.058X10" 4 

• C 56X10 -4 

.COXIO -4 

.063X10 -4 

.061X10" 4 

.063X10" 4 

.C65X10- 4 

.07OX10' 4 


12 

• 1C 5 

.126 

.129 

.163 

.160 

.134 

.140 

.151 

.144 

.168 

.160 

.1(7 

.190 

.175 

.181 

.191 

.199 


13 

.146 

.158 

.182 

.206 

.215 

.163 

.195 

.200 

.210 

.160 

.234 

.2*4 

.258 

.234 

.251 

.255 

.274 


14 

.C 13 

.005 

.010 

.014 

.010 

.014 

.003 

.010 

.017 

.010 

• CC7 

.CIO 

.005 

.007 

.006 

.010 

.011 

= 

15 

• C 24 

.026 

.023 

.024 

.022 

.029 

.029 

.027 

.024 

.026 

• C27 

• C 24 

.026 

.027 

.024 

.036 

.022 


16 

»C37 

.034 

.035 

.037 

.033 

.040 

.047 

.043 

.040 

.042 

.043 

. C 42 

.040 

.042 

.043 

.046 

.043 

— 

17 

«C 4t 

.040 

.045 

.045 

.048 

.050 

.060 

.055 

.061 

.056 

• C57 

.Cf 7 

.057 

.058 

.055 

.063 

.060 


18 

.C 72 

.070 

.070 

»C 76 

.074 

.080 

.088 

.083 

• Q B7 

.086 

• C 86 

.051 

.092 

.093 

.093 

.093 

.095 

= 

19 

• G <i l 

.099 

.103 

.103 

.113 

.104 

.123 

.123 

.115 

.123 

.129 

.135 

.157 

.141 

.143 

.144 

.153 

— 

20 

.163 

.179 

.194 

.219 

.227 

.198 

.220 

.231 

.213 

.243 

• 256 

.2(3 

.277 

.274 

.281 

.297 

.308 


21 

.2 ce 

.237 

.276 

.298 

.327 

*242 

• 293 

.311 

.323 

.346 

.368 

.329 

.407 

.371 

.304 

.414 

.435 

— 

22 

.2 31 

.280 

• 329 

.364 

.405 

.268 

.338 

.359 

.403 

.423 

.443 

.413 

.513 

.439 

,446 

.499 

.531 

— 

23 

• C47 

.049 

.050 

.049 

.047 

.055 

.059 

.057 

.065 

.058 

. C63 

.059 

.058 

.065 

• C6 7 

.065 

.066 


24 

*C!S 

.053 

.051 

.051 

.053 

• 062 

• 065 

.069 

.070 

.065 

.067 

.0(6 

.056 

.069 

.070 

.076 

.069 


25 

• C 14 

.073 

■ C70 

.069 

.068 

.082 

.091 

.087 

.079 

.087 

• C 65 

.063 

.083 

.093 

.093 

.095 

• 093 


26 

.0 « 

.089 

.089 

.093 

.094 

.105 

.112 

.116 

.122 

.11* 

.116 

.116 

.114 

.124 

.126 

.124 

.130 

— 

27 

.139 

.145 

.147 

.155 

.157 

.170 

.168 

.192 

.200 

.194 

.198 

.159 

• 2CC 

.220 

.223 

.226 

.223 

_ 

28 

• 2C6 

.2 12 

.231 

.239 

.252 

.192 

.294 

.289 

.303 

.309 

.318 

.329 

.329 

.428 

.353 

.360 

.375 

z 

29 

.2(6 

.319 

.357 

.393 

.423 

.354 

.420 

.438 

.460 

.480 

.501 

• 522 

.546 

.536 

.556 

1.096 

• 605 

s 

30 

.214 

.407 

.488 

.505 

.598 

.423 

.516 

.560 

.577 

.649 

• 6 97 

.722 

.757 

.697 

.732 

.796 

I .826 

- 

31 

.354 

.498 

• 608 

.707 

.760 

.512 

.611 

.679 

.710 

.787 

.641 

• 528 

.573 

.659 

• 920 

1.005 

! 1.063 


32 

• C77 

.077 

• C 71 

.079 

.071 

.091 

.098 

.094 

.088 

.089 

• C 88 

.066 

.093 

.106 

.102 

.103 

1 .101 


33 

.C*3 

.092 

.091 

.087 

.067 

.ire 

.120 

.116 

.122 

.113 

• 115 

.102 

.106 

.131 

.131 

.125. 

1 .125 

" 

34 

.146 

.144 

.139 

.147 

.147 

.095 

.192 

.108 

.134 

.123 

.186 

.1(6 

.112 

.212 

.214 

.212 

.216 

B 

35 

.223 

.226 

.228 

• 234 

.238 

.278 

.305 

.311 

.276 

.315 

.319 

.217 

.322 

.361 

.372 

.371 

.369 

- 

36 

.313 

.326 

.346 

.355 

.369 

• 400 

.453 

.457 

.451 

.480 

.491 

.504 

.505 

.558 

.563 

.574 

.587 


37 

1 .427 

.471 

.521 

■ .558 

' .597 

• 566 

1 .653 

.679 

.690 

.731 

.761 

i .165 

.803 

1 .839 

.860 

.902 

! .932 


r 38 

.521 

.617 

.721 

.608 

1 .878 

.699 

I .932 

.898 

.936 

1.006 

1.064 

1.122 

1.172 

[1.148 

1.193 

1.270 

1.337 


39 

.577 

.743 

.928 

1.080 

;1.209 

.817 

j .990 

1.096 

1.179 

1.294 

1.371 

1.4*7 

1.585 

1.456 

1.551 

1.670 

1.770 


40 

.617 

.817 

1.080 

11.294 

1.454 

• 871 

1.084 

1.234 

1.099 

1.507 

1.633 

1.756 

1.882 

1.663 

1.803 

1.960 

2.098 



.320 

.319 

.319 

1 .317 

1 .316 

.405 

1 .455 

.451 

.447 

.449 

.456 

.453 

• 444 

.533 

.531 

.528 

• 526 


1 42 

.3(1 

• 381 

• 388 

1 .393 

1 *397 

.493 

| .552 

.553 

.543 

.556 

.569 

.563 

.554 

1 .654 

.656 

• 664 

.666 


i 43 

.4 62 

• 468 

' .489 

.496 

.511 1 

.613 

.690 

.707 

.683 

.717 

.714 

.128 

.741 

I .838 

.844 

.853 

.871 


44 

.548 

.587 

, .630 

1 .658 

.688 

[ .764 

.876 

.898 

.871 

.931 

.952 

• 573 

• 986 

1.102 

1.119 

1.147 

1.177 


45 

.620 

.694 

; .776 

1 .836 

.891 

.895 

1.027 

1.075 

1.109 

1.159 

1.200 

1.252 

1.282 

.1*361 

1.408 

1.469 

1.521 


46 

.645 

• 636 

.572 

1.066 

1.187 

1.029 

1.216 

1.316 

1.347 

1.473 

1.550 

1.642 

1.703 

1.724 

1.816 

1.921 

2.009 


47 

• 7 59 

• 960 

1.208 

1.410 

1.592 

1.155 

1.430 

1.569 

1.657 

1.858 

1.978 

1.553 

2.244 

2.181 

2.325 

2.482 

2.623 


40 

.ece 

1.064 

'1.404 

1.786 

2.128 

1.262 

1.506 

1.810 

1.983 

2.239 

2.414 

2.6*6 

2.852 

12.613 

'2.831 

3.074 

3.299 


49 

.667 

1.173 

1.594 

2.132 

2.803 

1.363 

1.753 

2.018 

2.514 

2.579 

2.829 

3.1C5 

3.359 

,3.027 

3.288 

3.612 

3.902 


50 

.814 

1.00 8 

1.143 

1.245 

1.347 

1.440 

1.618 

1.739 

1.759 

1.832 

1.550 

2 .025 

2.090 

2.330 

2.360 

2.483 

2.724 


51 

.859 

1.093 

1.302 

1.424 

1.570 

1.528 

1.774 

1.900 

1.931 

2.133 

2.253 

<2.356 

2.428 

2.643 

2.765 

2.912 

3.026 


52 

1 .454 

1.159 

1.414 

1.636 

1.851 

1.577 

1.928 

2.136 

2.225 

2.451 

2.600 

2.757 

2.905 

3.040 

3.216 

3.415 

3.617 


53 

• 462 

1.236 

1.557 

1.844 

2.092 

1.637 

2.011 

2.271 

2.471 

2.722 

2.514 

3.150 

3.332 

3.450 

3.649 

3.094 

4.096 


54 

.554 

1.292 

1.602 

,1.964 

2.262 

1.679 

2.067 

2.351 

2.590 

2.901 

3.135 

3.425 

3.633 

1.713 

3.961 

4.279 

4.528 

B 

55 

1.058 

1.393 

1.824 

12.230 

2.624 

1.806 

2.253 

2.608 

2.914 

3.241 

3.612 

3.535 

4.236 

4.201 

4.534 

4.072 

5.293 

“ 

56 

1.1C9 

1.460 

1.941 

2.412 

2.846 

1.881 

2.392 

2.766 

3.100 

3.543 

3.905 

4.265 

4.677 

4.566 

5.007 

5.485 

5.908 


57 

1.139 

.955 

2.060 

2.588 

3.059 

1.930 

2.496 

2.903 

3.258 

3.790 

4.172 

4.7C2 

5.211 

4.882 

5.427 

6.029 

6.578 


58 

1.215 

1.597 

1 1**578 

2.379 

2.751 

2.161 

2.719 

3.081 

3.381 

3.849 

1 4.133 

4.6(7 

4.927 

5.716 

6.142 

6.526 

6.891 

■ 

59 

1.265 

1.676 

2.150 

2.621 

3.042 

2.260 

2.870 

3.282 

3.575 

4.186 

4.562 

5.113 

5.498 

i6. 165 

7.019 

7.615 

8.319 

1 

60 

1.268 

1.729 

j2.275 

2.800 

3.287 

2.357 

2.941 

3.446 

3.851 

4.470 

4.671 

5.462 

6.C72 

6.528 

7.566 

8.750 

9.833 


61 

1.346 

1.806 

2.407 

3.012 

3.532 

2.391 

3.070 

3.615 

4.121 

4.769 

5.270 

5.563 

6.623 

6.932 

8.200 

9.797 

11.652 
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TABLE IV 

STRUCTURAL INFLUENCE -COEFFICIENT MEASUREMENTS - Continued 
(b) Matrix [b] - Continued 


Load 

Panel deflection per unit load, meter s/newt on, measured at position: 


point 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

II 

.05flXl(T 4 

.061XKT 4 

•064X10- 4 

.071x10-4 

.026X10' 4 

• C30XKT 4 

.025X10-4 

.042x10-4 

.051X10-4 

.064X10 -4 

.070X10-4 

.077X10- 4 

.075X10-4 

•039X10" 4 

.051X10-4 

•049X10 -4 

• Q 66X1 O'* 4 

12 

• 078 

.099 

• 109 

.126 

.028 

.034 

.034 

.063 

.074 

• 100 

.119 

.122 

.128 

• 053 

• 064 

• 073 

,007 

13 

• 116 

.150 

• 181 

.212 

.029 

• 048 

• 067 

.083 

.106 

.161 

.195 

.216 

.225 

• 076 

• 094 

• 111 

.136 

14 

.002 

.002 

• 003 

.004 

-.001 

• 002 

• CC3 

.005 

-.003 

.002 

.001 

-.010 

.005 

-.003 

-.005 

-.007. 

-.006 

15 

• 009 

• 008 

.007 

.007 

.003 

.003 

• CC6 

• 003 

.003 

•009 

.006 

.010 

.004 

• 003 

•005 

.007 

.006 

16 

.012 

• 013 

.013 

• 014 

.006 

.005 

• Cll 

• 009 

.002 

•Oil 

• 014 

-.014 

.013 

• 005 

.010 

.008 

• 002 

17 

• 023 

.025 

.023 

• 026 

.013 

.019 

.015 

• 017 

.015 

•026 

.026 

.017 

.010 

• 016 

.017 

.018 

•C09 

18 

.030 

• 034 

• 034 

• 038 

.017 

.016 

•021 

• 024 

.019 

• 033 

• 036 

.029 

.035 

• 018 

.015 

• 023 

.019 

19 

.042 

.052 

.056 

.061 

.018 

.021 

•on 

.027 

.029 

.042 

.053 

.047 

.057 

• 019 

.021 

.026 

• 035 

20 

.059 

.077 

• 085 

.097 

.019 

.024 

.022 

.030 

.047 

.069 

.085 

.078 

.094 

.029 

.031 

• 043 

• 044 

21 

• 072 

• 095 

.107 

.124 

.024 

.033 

.043 

.04 3 

.047 

.005 

.107 

.117 

.125 

• 031 

.038 

• 050 

.059 

22 

• 062 

• 006 

.090 

.115 

.011 

• C19 

.017 

.030 

.030 

.076 

• 096 

.099 

.112 

.017 

.023 

• 033 

• C46 

23 

• 005 

.005 

.006 

.007 

-.004 

-.003 

• C03 

-.002 

• 008 

.003 

• 005 

-.007 

.008 

-.002 

-.004 

-.005 

-•C03 

24 

• 007 

• 009 

• 010 

• 010 

-.002 

-.001 

-.003 

• 003 

• 0C9 

•003 

.010 

.009 

.009 

-.003 

-.003 

-.002 

.012 

25 

• 010 

.011 

.011 

.012 

.001 

.003 

-.015 

•00 5 

.001 

.005 

.005 

.006 

.005 

-.003 

-•005 

-•002 

-.002 

26 

• 014 

.017 

• 010 

.020 

.003 

.002 

• CC2 

• 004 

0.000 

•Oil 

.010 

.013 

.015 

-.008 

-•005 

-.006 

-.005 

27 

• 019 

• 026 

• 029 

.033 

.004 

.008 

-.001 

.003 

.003 

.010 

.019 

.024 
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TABLE V 

LONGITUDINAL AERODYNAMIC CHARACTERISTICS - INFLUENCE -COEFFICIENT METHOD 
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Figure 4.- Flexible model prior to wind-tunnel test. 
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Figure 6.- Model control points for vibration test, 
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Figure 7.- Continued. 
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Figure 8.- Test setup for measuring structural influence coefficients. 
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Figure 9.- Control -point location and paneling scheme for influence -coefficient analysis. 
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Figure 11.- Effect of dynamic pressure on measured and calculated 
flexible -to -rigid Ct and aerodynamic -center movement. 
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Figure 11.- Continued. 
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Figure 11.- Continued. 
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(a) M = 0.6. 

Figure 14.- Effect of dynamic pressure on lateral -directional stability derivatives. 
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(c) M = 1.2. 
Figure 14.- Continued. 
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Figure 14.- Continued. 
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